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THE  EFFECT  OF  FLUIDS  ON  WAVEGUIDES  BELOW  CUTOFF 
PENETRATIONS  AS  RELATED  TO  ELECTROMAGNETIC 
SHIELDING  EFFECTIVENESS 


1  INTRODUCTION 


Background 

It  is  common  practice  in  electromagnetic  pulse  (EMP)  and  electromagnetic 
interference  (EMI)  hardened  facilities  that  house  sensitive  electronic 
equipment  to  enclose  part  or  ail  of  the  facility  in  a  metallic  electromagnetic 
shield.  This  shield  must  have  openings  for  personnel  and  equipment  access, 
electric  power  entry,  communications  and  control  cooling,  heating, 
ventilating,  air-conditioning  (HVAC),  and  entry  of  various  fluids  such  as 
water  and  refrigerants.  Technology  for  the  design  of  these  penetrations  is, 
in  general,  highly  advanced  and  is  discussed  in  numerous  technical  manuals  and 
design  handbooks.1 

When  it  is  required  to  bring  fluids  through  a  shield,  it  must  be  done  in 
a  manner  which  does  not  compromise  the  shield.  One  approach  is  to  use  a 
completely  closed  metallic  piping  or  conduit  system  into  which  compromising 
electromagnetic  energy  cannot  enter.  The  metallic  pipe  or  conduit  typically 
is  welded  or  otherwise  electrically  connected  around  its  periphery  to  the 
shield  at  the  point  of  entry. 

Another  approach  is  to  use  a  short  length  of  metallic  pipe  to  penetrate 
the  shield.  This  length  of  pipe  acts  as  a  waveguide  below  cutoff,  blocking 
unwanted  electromagnetic  energy.  (A  waveguide  passes  only  those  frequencies 
above  a  particular  "cutoff"  frequency.)  This  technique  is  used  to  isolate  the 
pipe  ground  from  the  shield  ground  without  compromising  the  shield.  It  also 
permits  the  connection  of  nonmetallic  pipes  to  the  shield  without  a  loss  in 
shielding. 

A  honeycomb  grating  is  an  extension  of  the  waveguide  below  cutoff 
concept.  A  honeycomb  grating  is  made  of  many  tubular  cells  with  hexagonal 
cross  sections.  Each  cell  acts  as  a  waveguide  below  cutoff. 

Typically,  when  designing  shield  penetrations,  the  pipes  are  assumed  to 
be  air-filled.  In  practice,  however,  it  is  often  necessary  to  bring  fluids 
other  than  air  through  the  pipes.  Examples  are  tapwater,  distilled  water, 
fuel,  refrigerants,  fire  extinguishing  fluids,  and  freeze  prevention  fluids. 
These  fluids  may  have  relative  dielectric  constants  and  dissipation  factors 


DASA  EMP  Handbook,  DASA  2114-1  (Defense  Atomic  Support  Agency  [DASA]  Infor¬ 
mation  and  Analysis  Center,  September  1968);  DNA  EMP  Awareness  Course  Notes, 
Third  Ed,  DNA  2772T  (Defense  Nuclear  Agency,  October  1977);  EMP  Engineering 
Practices  Handbook,  NATO  file  No.  1460-2  (October  1977);  Nuclear  Electro¬ 
magnetic  Pulse  (NEMP)  Protection,  TM  5-855-5  (Headquarters,  Department  of  the 
Army,  February  1974). 
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much  different  from  air.  The  attenuation  of  the  electromagnetic  energy 
propagating  within  the  fLuid-tilled  pipe  will  be  different  from  that  of  the 
air-filled  pipe. 


In  acceptance  testing  of  EMI  and  EMP  hardened  facilities,  the  pipe 
penetrations  usually  are  checked  using  the  techniques  defined  in  MIL-STD  285 
or  IEEE  299.2  These  methods  use  horn  antennas  or  open  waveguide  antennas  to 
measure  the  shielding  effectiveness.  Measurement  is  relatively 
straightforward  for  a  simple  shield  but  becomes  more  complicated  when  the 
shield  has  penetrations.  Some  antenna  positions  required  by  IEEE  299  cannot 
be  achieved  because  the  external  piping  connected  to  the  shield  prevents 
positioning  of  the  test  antennas  on  the  pipe  axis.  The  same  is  true  for  some 
honeycomb  gratings. 

Clearly,  fluid-filled  pipes  penetrating  a  shielded  enclosure  have  special 
design  requirements.  Designers  of  EMI  or  EMP  shields  therefore  need  specific 
data  to  assist  in  designing  fluid  penetrations. 


Objective 

The  objectives  of  this  study  were  to: 

1.  Perform  a  mathematical  analysis  of  waveguides  and  honeycomb  filters 
to  determine  attenuation  versus  frequency. 

2.  Provide  design  curves  showing  cutoff  frequency  and  attenuation  versus 
frequency  for  commonly  used  pipe  sizes  and  fluids. 

3.  Verify  the  above  solutions  empirically. 


Approach 

A  literature  search  was  performed  to  obtain  information  on  related  work 
covering  energy  capture  ratios  for  open  waveguides,  measurement  of 
electromagentic  parameters  of  fluids,  and  performance  of  honeycomb  filters. 
Computer  programs  were  developed  to  calculate  the  waveguide's  performance  when 
fluid-filled.  Laboratory  testing  was  performed  on  representative  waveguides 
to  verify  analytical  methods.  Analysis  was  performed  to  calculate  the 
shielding  effectiveness  of  honeycomb  filters. 


Mode  of  Technology  Transfer 

It  recommended  that  the  results  of  this  study  be  incorporated  in  the 
Defense  Nuclear  Agency's  EMP  Engineering  Practices  Handbook  for  Facilities 


2Method  of  Attenuation  Measurements  for  Enclosures,  Electro-Magnetic 
Shielding,  for  Electronic  Test  Purposes,  MIL-STD  285  (Department  of  Defense, 
July  1956);  Measurement  of  Shielding  Effectiveness  of  High-Performance 


which  is  in  preparation.  it  is  also  recommended  that  these  results  be 
included  in  Army  Technical  Manual  (TM)  5-855-5,  Nuclear  Electromagnetic  Pulse 
( NEMP)  Protection . 


2  THEORY  AND  PROCEDURE 


Introduction 

This  chapter  describes  the  anaLysis  of  a  dielectric  loaded  waveguide. 
Specifically,  the  modes  excited  in  a  waveguide  by  an  incident  plane  wave  and 
the  accuracy  of  shielding  tests  such  as  IEEE  299  and  MIL-STD-285  are 
addressed.  An  experimental  verification  for  the  theoretical  analysis  of  the 
loaded  waveguide  is  included. 

A  waveguide  is  a  hollow  metallic  pipe  (usually  circular  or  rectangular) 
through  which  electromagnetic  energy  propagates.  The  most  important  waveguide 
parameter  is  its  cutoff  frequency,  that  is,  the  frequency  at  which  the 
waveguide  starts  dissipating  the  electromagnetic  energy  instead  of  propagating 
it.  Frequencies  above  cutoff  are  passed  through  the  waveguide  with  no 
attenuation.  Frequencies  beLow  cutoff  are  heavily  attenuated  within  the 
waveguide. 

The  waveguide  cutoff  frequency  is  a  function  of  the  waveguide  cross 
sectional  area  and  the  contents  of  the  waveguide.  The  larger  the  cross 
sectional  area,  the  lower  the  cutoff  frequency.  Generally,  any  fluid  in  the 
waveguide  besides  air  will  have  a  lower  cutoff  frequency  than  the  air-filled 
waveguide. 

A  waveguide  below  cutoff  penetration  of  a  shield  consists  of  a  length  of 
metal  pipe  passing  through  the  shield.  Where  the  pipe  penetrates  the  shield, 
the  pipe  periphery  is  electrically  connected  to  the  shield. 

Waveguide  below  cutoff  penetrations  are  used  in  plastic  pipe  systems  to 
block  unwanted  electromagnetic  energy.  They  are  also  used  to  eliminate  ground 
loops . 

In  most  shielded  facilities,  numerous  piping  and  conduit  penetrations 
will  be  required,  possibly  at  widely  spaced  points.  The  piping  and  conduit 
outside  the  shield  act  as  electromagnetic  energy  collectors.  Many  of  these 
pipes  and  conduits  pass  underground  and  thus  are  earth-grounded.  The  use  of 
many  earth-grounded  piping  systems  with  shield  penetrations  creates 
electromagnetic  energy  collecting  ground  loops.  These  ground  loops  are  not 
consistent  with  optimized  grounding  system  philosophy,  so  it  may  be  necessary 
to  electrically  disconnect  the  piping  system  from  the  shield.  This  is 
commonly  done  by  inserting  nonmetallic  pipe  sections  near  the  shield,  then 
penetrating  the  shield  with  a  waveguide  below  cutoff  piping  section. 


* 


t 


Circular  Waveguide  Theor\ 


Figure  1  shows  the  basic  geometry  of  a  pipe  penetrating  a  shielded 
enclosure.  The  pipe  has  inside  diameter  d  and  length  L.  It  contains  a  fluid 
with  relative  dielectric  constant  e  and  .oss  tangent  tan  6.  A  plane  wave  is 
incident  from  the  left,  propagating  along  the  z-axis. 
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A  rigorous  analysis  of  the  electromagnetic  energy  coupled  into  the  pipe 
and  then  into  the  enclosure  is  difficult.  Because  the  relative  dielectric 
constants,  loss  tangents,  frequencies,  and  pipe  dimensions  vary  widely, 
several  types  of  analysis  are  required.  Most  importantly,  a  rigorous 
multimodal  examination  of  the  problem  would  require  precise  knowledge  of  the 
excitation  mechanism. 

The  forms  of  excitation  found  in  practice  are  likely  to  vary  widely  and 
would  most  certainly  differ  from  the  basic  geometry  just  described.  Because 
of  this  difference,  an  approximate  transverse  electric  and  magnetic  fields 
(TEM)-waveguide  analysis  is  used  to  obtain  the  shielding  effectiveness  of  the 
pipe  alone.  Inclusion  of  a  particular  excitation  mechanism  would  increase 
shielding  effectiveness  because  the  coupling  of  the  incident  wave  to  the  pipe 
is  imperfect. 

Appendix  A  presents  the  basic  equations  for  the  analysis  of  electro¬ 
magnetic  wave  propagation  through  a  circular  waveguide  (pipe). 


Modes  Excited  in  a  Circular  Waveguide 

Typically,  most  of  the  energy  in  a  waveguide  is  carried  in  the  first  few 
propagating  modes.  This  was  observed  during  the  tests  conducted  for  this 
report . 

Figure  2  is  a  plot  of  energy  propagation  through  a  1.5-in.  (3.81-cm) 
inside  diameter  (i.d.)  pipe  as  a  function  of  frequency.  The  pipe  was  excited 
by  a  plane  wave  as  shown  in  Figure  1.  The  vertical  scale  represents  the 
attenuation  resulting  from  the  pipe  relative  to  free  space. 
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Figure  2.  Transmission  through  an  air-filled  1.5-in.  (3.8-cm) 
i.d.  pipe. 


P 


For  an  air-filled  pipe,  the  two  lowest  modes  are  the  TE,,  and  TMq,  modes 
with  cutoff  frequencies  of  4.7  and  6.14  GHz,  respectively.  Figure  2  clearly 
shows  that  the  TEj,  must  be  present.  Also,  there  is  no  significant  increase 
in  the  power  transferred  above  the  TMq^  cutoff  frequency,  indicating  that  most 
of  the  energy  is  coupled  into  Che  TEj^  mode. 


Computer  Model  of  a  Circular  Waveguide 

Appendix  C  contains  listings  and  descriptions  of  two  FORTRAN  V  programs,  ~ 

Pipe  1A  and  Pipe  IB  (Figures  Cl  and  C2),  that  compute  the  theoretical 
attenuation  resulting  from  a  length  of  pipe.  Both  programs  are  designed  for 
use  on  an  interactive  computer  system.  The  frequency  and  corresponding 
attenuation  are  written  into  a  named  file  for  storage.  The  user  is  questioned 
for  the  problem  parameters.  ^ 

Program  Pipe  1A  computes  the  attenuation  over  a  band  of  frequencies.  The 
relative  dielectric  constant  and  loss  tangent  are  assumed  to  be  constant  over 
this  band.  The  variables  to  be  entered  are:  pipe  inside  radius,  pipe  length, 
relative  dielectric  constant,  loss  tangent,  minimum  frequency  in  band,  maximum 
frequency  in  band,  the  frequency  increment  to  be  used  within  the  band,  and  the  0 

name  of  the  output  file  in  which  the  results  are  to  be  stored. 

Program  Pipe  IB  computes  the  attenuation  for  a  particular  frequency.  The 
program  is  cyclical  so  the  attenuation  at  any  number  of  frequencies  may  be 
computed.  However,  the  pipe  inside  radius  and  pipe  length  are  entered  only 
once.  The  parameters  to  be  entered  are:  pipe  inside  radius,  pipe  length,  the  9 

name  of  the  output  file  in  which  the  results  are  to  be  stored,  frequency  1, 
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relative  dielectric  constant  1,  loss  tangent  1,  frequency  2,  relative 
dielectric  constant  2,  loss  tangent  2,  frequency  3,  .... 


Experimental  Verification 

A  crude  verification  of  the  theoretical  waveguide  attenuation  was 
performed.  The  experimental  set-up  is  shown  in  Figure  3.  The  system  consists 
of  a  6-in.  (15.25-cm)  pipe  penetrating  a  shielded  enclosure.  The  pipe  inside 
diameter  is  1.48  in.  (3.76  cm).  The  ends  of  the  pipe  were  covered  with 
plastic  wrap  so  the  pipe  could  be  filled  with  fluid.  The  receiver  and  its 
horn  antenna  were  placed  inside  the  shielded  enclosure.  The  signal  generator 
and  transmitting  horn  (identical  to  the  receiving  horn)  were  placed  outside 
the  enclosure.  The  distance  between  each  horn  and  the  nearest  end  of  the  pipe 
was  12  in.  (30.48  cm),  resulting  in  an  overall  horn  separation  of  30  in. 

(76.20  cm)  (12  +  6  +  12).  Free-space  reference  measurements  were  taken  with 
both  horns  outside  the  enclosure  spaced  30  in.  (76.20  cm)  apart. 

Measurements  were  taken  between  2.5  and  7.5  GHz  with  a  0.5-GHz  spacing. 

An  additional  measurement  was  made  at  7.6  GHz.  This  narrow  test  band  was  a 
result  of  equipment  limitations.  (The  horns  could  not  be  used  below  2.5  GHz 
and  the  transmitter  could  not  be  used  above  7.6  GHz.) 

Incident  wave  power  is  lost  not  only  in  the  waveguide  but  also  at  the 
free  space-pipe  interfaces.  The  power  is  lost  at  these  interfaces  because  of 
coupling.  Therefore,  to  verify  the  waveguide  attenuation,  the  coupling  loss 
must  also  be  determined. 

It  is  best  to  measure  the  coupling  loss  as  it  represents  the  loss  due  to 
a  particular  waveguide  excitation  mechanism.  Also,  since  excitation 
mechanisms  vary  widely  depending  on  geometry,  nearby  metallic  objects,  and 
wave  type,  it  is  more  practical  to  measure  coupling  loss. 

The  measurement  of  coupling  loss  is  simple  for  frequencies  above 
cutoff.  The  wave  will  not  be  attenuated  within  an  air-filled  pipe. 

Therefore,  the  difference  between  the  free  space  reference  measurement  and  the 
measurement  for  the  air-filled  pipe  will  be  the  coupling  loss. 

For  frequencies  below  cutoff,  the  measurement  of  coupling  loss  is 
dependent  on  the  wave  attenuation  within  the  pipe.  For  this  reason,  a 
theoretical  method  was  used  to  determine  the  coupling  loss.  Appendix  B 
describes  the  calculation  of  coupling  loss  using  antenna  directivities.  Table 
1  lists  the  computed  and  measured  coupling  losses  for  the  experimental  set-up 
of  Figure  3. 

Four  fluids  covering  a  wide  range  of  relative  dielectric  constants  were 
used  in  the  test:  air  (e  =  1),  carbon  tetrachloride  (e  »  2),  meth¬ 
anol  (e  »  24),  and  distilled  water  (e  ~  63).  Tables  5  through  4  list 
resultsrof  the  test.  For  methanol  andrdistilled  water,  the  frequency  band  was 
limited  because  of  high  waveguide  attenuation  (-90  dBm  was  the  minimum 
receiver  sensitivity). 
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Table  2 


Transmission  Loss  Through  1.5-In.  i.d. 

With  Carbon  Tetrachloride 

Pipe  Filled 

Frequency 

(GHz) 

Free-space 

reference 

(dB) 

Theoretical 
attenuation 
in  pipe 
(dB) 

Coupling 

loss 

(dB) 

Computed 

measurement 

Actual 

measurement 

2.5 

-5 

81 

24 

-110 

-79 

3.0 

13 

44 

22 

-53 

-37 

3.5 

11 

0 

21 

-10 

-  8 

4.0 

7 

0 

20 

-13 

-  8 

4.5 

7 

0 

19 

-12 

-  8 

5.0 

7 

0 

18 

-11 

-11 

5.5 

2 

0 

14 

-12 

-12 

6.0 

12 

0 

24 

-12 

-13 

6.5 

-7 

0 

11 

-18 

-19 

7.0 

-19 

0 

13 

-32 

-33 

7.5 

-15 

0 

24 

-39 

-38 

7.6 

-11 

0 

32 

-43 

-41 

Table  3 

Transmission  Loss  Through  1.! 

5-In,  i.d.  Pipe 

Filled  With  Methanol 

Reference 

Frequency 

free 

Coupling 

Theoretical 

Computed 

Actual 

(GHz) 

space 

loss 

attenuation 

measurement 

measurement 

4.0 

7 

20 

65 

-78 

-83 

4.5 

7 

19 

67 

-79 

-87 

5.0 

7 

18 

65 

-76 

-81 

Table  4 

Transmission  Loss 

Through  1.5- In.  i.d.  Pipe  Filled  With  Distilled  Water 

Reference 

Frequency 

free 

Coupl ing 

Theoretical 

Computed 

Actual 

(GHz) 

space 

loss 

attenuation 

measurement 

measurement 

2.5 

-5 

24 

47 

-76 

-74 

3.0 

13 

22 

58 

-67 

-72 

3.5 

11 

21 

68 

-78 

-79 

4.0 

7 

20 

67 

-80 

-77 

4.5 

7 

19 

66 

-78 

-84 

5.0 

7 

18 

59 

-70 

-82 

5.5 

2 

14 

54 

-66 

-85 

6.0 

12 

24 

53 

-65 

-79 

6.5 

-7 

11 

52 

-70 

-84 
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The  theoretical  attenuation  in  the  pipe  was  found  using  program  Pipe 
IB.  Because  of  the  lack  of  published  data  on  relative  dielectric  constants 
and  loss  tangents,  these  parameters  were  measured  using  a  dielectric  slab 
technique.  Appendix  D  describes  this  technique  and  lists  the  values  of  the 
relative  dielectric  constants  and  loss  tangents  used  in  computing  the 
waveguide  attenuation. 

The  predicted  measurement  is  equal  to  the  free-space  reference  minus  the 
sum  of  the  coupling  loss  and  theoretical  waveguide  attenuation  (i.e.,  column  5 
=  column  2  -  (column  3  +  column  4j  in  Tables  2  through  4).  The  Last  coLumn 
contains  the  actual  measurement. 

For  carbon  tetrachloride,  the  agreement  between  predicted  and  actual 
measurements  is  fairly  good  above  cutoff.  Below  cutoff,  the  error  is 
extremely  large.  No  reason  for  this  disagreement  is  known. 

For  methanol  and  distilled  water,  the  actual  measurement  is  consistently 
less  than  the  predicted  measurement  (except  in  two  cases).  This  can  be 
explained  by  the  reflection  and  refraction  of  the  electromagnetic  wave  at  the 
air-fluid  interface.  Electromagnetic  waves  undergo  reflection  and  refraction 
at  interfaces  where  the  relative  dielectric  constants  change.  The  reflection 
coefficient,  R,  is  given  by: 


R  = 


^1-^2 

/Ti  ♦  /r2 


(Eq  II 


where  is  the  relative  dielectric  constant  of  the  first  medium  and  z^  is  the 
relative  dielectric  constant  of  the  second  medium.  Therefore,  the  greater  the 
mismatch  of  relative  dielectric  constants,  the  greater  the  reflection  at  the 
interface.  This  means  there  will  be  an  additional  power  loss  due  to 
reflection  and  thus,  a  lower  value  of  transmitted  power. 


Shielding  Tests  for  Circular  Waveguides 

It  is  important  to  consider  the  accuracy  of  standard  shielded  enclosure 
test  procedures  in  measuring  the  shielding  effectiveness  of  enclosures  with 
pipe  penetrations.  The  cutoff  frequency  of  a  pipe  is  inversely  proportional 
to  the  pipe  radius  and  to  the  square  root  of  the  relative  dielectric  constant 
of  the  fluid  in  the  pipe  (i.e.,  as  r  or  /z~  increases,  the  cutoff  frequency 
decreases).  Neither  MIL-STD-285  nor  IEEE  $99  make  any  statement  about  the 
contents  of  the  penetrating  pipe  whiLe  it  is  being  tested.  Since  any  fluid 
besides  air  in  the  pipe  will  Lower  the  cutoff  frequency,  it  is  important  to 
have  the  pipe  fiLLed  as  it  would  be  under  normal  use.  Since  the  greatest 
amount  of  electromagnetic  energy  is  coupled  into  the  pipe  when  the  plane  wave 
is  axially  incident,  the  shielding  effectiveness  measurements  should  be  taken 
with  the  transmitting  and  receiving  antennas  positioned  along  the  pipe  axis. 


3  RESULTS 


This  chapter  presents  a  series  of  graphs  produced  from  FORTRAN  Program 
Pipe  IB  to  be  used  in  designing  pipe  penetrations  for  shielded  enclosures. 
Before  presenting  these  graphs,  a  few  comments  will  be  made  about  the  effect 
of  material  parameters  on  waveguide  attenuation  and  cutoff  frequency. 


Effect  of  Material  Parameters  on  Attenuation 


The  loss  tangent  of  the  dielectric  affects  attenuation  whereas  the 
relative  dielectric  constant  affects  cutoff  frequency.  By  definition,  the 
cutoff  frequ  ,ncy  affects  the  attenuation.  Figure  4  gives  the  attenuation  of 
the  lowest  order  mode  (TEj,)  of  a  6-in.  (15.24-cm)  length  of  1.475-in.  (3.746- 
cm)  i.d.  pipe  for  distilled  water  at  various  loss  tangents.  The  loss  tangent 
and  the  relative  dielectric  constant  are  held  constant  throughout  the  , 
frequency  band.  The  value  of  the  relative  dielectric  constant  is  68.11.  The 
loss  tangent  is  stepped  from  0.0  to  0.50  in  0.05  increments.  The  major  effect 
of  the  loss  tangent  is  to  increase  the  attenuation.  Note  that  a  small  change 
in  the  loss  tangent  can  change  the  attenuation  considerably,  especially  at 
frequencies  well  above  cutoff. 

Curves  showing  the  change  in  cutoff  frequency  of  the  lowest  order  mode 
with  relative  dielectric  constant  are  given  in  Figure  5.  The  curves  represent 
pipe  i.d.  of  1.0  in.  (2.54  cm),  1.5  in.  (3.81  cm),  2.0  in.  (5.08  cm),  2.5  in. 
(6.35  cm),  4.0  in.  (10.16  cm),  and  6.0  in.  (15.24  cm).  The  relative 
dielectric  constants  range  from  1  to  100.  Note  the  steepness  of  each  curve 
for  low  relative  dielectric  constants.  This  is  caused  by  the  inverse  square 
root  dependence  of  cutoff  frequency  on  the  relative  dielectric  constant. 


Design  Curves 

Figures  6  through  14  give  the  attenuation  of  the  lowest  order  mode  (TE^) 
in  a  6-in.  (15.24-cm)  length  of  pipe  for  various  fluids  and  pipe  i.d. 
values.  Each  figure  consists  of  six  curves  corresponding  to  i.d.  values  of 
1.0  in.  (2.54  cm),  1.5  in.  (3.81  cm),  2.0  in.  (5.08),  2.5  in.  (6.35  cm),  4.0 
in.  (10.16  cm),  and  6.0  in.  (15.24  cm).  The  relative  dielectric  constants  and 
loss  tangents  were  taken  from  von  HippelJ  and  are  summarized  in  Table  5.  For 
intermediate  frequencies,  a  linear  interpolation  was  performed  on  the  given 
values.  Because  attenuation  is  directly  proportional  to  length,  to  find  the 
attenuation  from  a  pipe  of  length  L,  the  attenuation  should  be  multiplied  by 
L/6  (L  in  in.). 

A  few  liquids  have  large  changes  in  relative  dielectric  constant  and  loss 
tangent  over  the  100-MHz  to  10-GHz  range  (e.g.,  methanol,  ethylene  glycol,  and 
water).  This  might  present  problems  in  determining  the  optimal  pipe  size, 
especially  since  a  small  change  in  loss  tangent  can  result  in  a  large  change 
in  attenuation  (see  Figure  4). 


JA.  R.  von  Hippel,  Dielectric  Materials  and  Applications  (Technology  Press  of 
the  Massachusetts  Institute  for  Technology,  1954). 
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'igure  A.  Theoretical  attenuation  of  the  TE^  mode  for  a  6-in. 

(15.2  cm)  length  of  1.5-in.  (3.8-cm)  i.d.  pipe  filled 
with  distilled  water  for  various  loss  tangents. 


Figure  5.  Cutoff  frequency  versus  relative  dielectric  constant  for 
various  pipe  diameters. 
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Figure  8.  Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  filled 
with  aviation  gasoline  (91  octane)  (e  *>  2)  for  various 
inside  diameters. 
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Figure  9.  Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  filled 

with  iet  fuel  (JP-3)  (e  a  2)  for  various  inside  diameters. 
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Figure  10.  Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  filled 
with  carbon  tetrachloride  (e^  =»  2)  for  various  inside 
diameters. 
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Figure  11.  Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  filled 
with  cable  oil  (e  =  2)  for  various  inside  diameters. 
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gure  12.  Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  filled 

with  methanol  (e  «  24)  for  various  inside  diameters, 
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Figure  13.  Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  filled 

with  ethylene  glycol  (e^  »  25)  for  various  inside  diameters. 
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Figure  14. 
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Theoretical  attenuation  of  a  6-in.  (15.25-cm)  pipe  "filled 
with  distilled  water  (e  »  63)  for  various  inside  diameters, 


A  worst-case  analysis  could  be  used  in  these  cases.  For  example,  a 
worst-case  condition  could  assume  a  loss  tangent  of  zero  and  a  large  relative 
dielectric  constant  (i.e.,  one  greater  than  or  equal  to  the  maximum  possible 
relative  dielectric  constant).  Hdwever,  this  could  lead  to  problems  in 
practice.  In  designing  a  pipe  penetration  to  carry  distilled  water,  for 
example,  assuming  a  relative  dielectric  constant  of  80,  the  maximum  cutoff 
frequency  (from  Figure  5)  is  approximately  750  MHz  for  a  1.0-in.  (2.54-cm) 
i.d.  pipe.  With  a  loss  tangent  of  zero,  little  attenuation  will  be  expected 
above  cutoff.  To  raise  the  cutoff  frequency,  either  a  smaller  diameter  pipe 
is  needed  or  a  honeycomb  grating  can  be  placed  in  the  pipe  (Appendix  E). 

Either  method  will  cause  a  reduction  in  the  flow  rate  of  the  fluid  in  the 
pipe.  The  other  alternative  is  to  assume  some  loss  tangent;  then  even  at 
cutoff  some  attenuation  will  be  present.  For  the  values  of  Table  5,  Figure  14 
indicates  an  attenuation  of  approximately  10  dB  at  cutoff  for  the  1-in.  (2.54- 
cm)  i.d.  pipe.  Therefore,  to  obtain  70  dB  of  attenuation,  a  pipe  of  length  L, 
with  L  given  by 


70  dB  _  L 
10  dB  6  in. 


would  be  required.  In  this  case,  L  =  42  in.  —  a  large  length  for  a  pipe 
penetration. 

This  calculation  is  based  on  values  interpolated  from  from  the  data  in 
Table  5.  If  the  loss  tangent  and  relative  dielectric  constant  were  known  in 
greater  detail,  at  smaller  frequency  increments,  perhaps  a  more  practical 
design  could  be  achieved. 
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Relative  Dielectric  Constants  and  Loss  Tangents  for 
Various  Materials* 


Frequency  Relative  dielectric  constant  Loss  tangent 

(MHz)  (tf)  (tan  6) 


Air 


100 

1 

0 

1000 

1 

0 

10000 

1 

0 

Aviation  gasoline  -  100  Octane 

300 

1.94 

.00008 

3000 

1.92 

.0014 

Aviation  gasoline  -  91  Octane 

300 

1.95 

.00004 

3000 

1.94 

.0015 

Jet  fuel  -  JP-3 

300 

2.08 

.0007 

3000 

2.04 

.0055 

Carbon  tetrachloride 

100 

2.17 

<-0002 

300 

2.17 

<.0001 

3000 

2.17 

.0004 

10000 

2.17 

.0016 

Cable  oil  -  Type  5314  (CE) 

300 

2.24 

.0039 

3000 

2.22 

.0018 

10000 

2.22 

.0022 

Methanol  at  25°C 

100 

31.0 

.038 

300 

30.9 

.080 

3000 

23.9 

.640 

10000 

8.9 

.810 

Ethlyene  Glycol  at  25'C 

100 

41 

.045 

300 

39 

.160 

3000 

12 

1.000 

10000 

7 

.780 

Distilled  water  at  25*0 

300 

77.5 

.016 

3000 

76.7 

.157 

10000 

55.0 

.540 

•From  A.  R.  von  HippeL.  These  data  were  used  to  plot  the  curves  shown  in 
Figures  6  through  14. 


4  CONCLUSIONS 


This  report  has  addressed  the  electromagnetic  shielding  effectiveness  of 
waveguide  sections  (pipes)  filled  with  fluids.  An  analytical  method  for 
determining  the  coupling  between  horn  antennas  when  used  to  measure  shielding 
effectiveness  of  shields  with  pipe  penetrations  has  been  presented. 

Analytical  methods  including  computer  codes  for  calculating  attenuation  versus 
frequency  within  the  pipe  sections  have  been  developed  and  verified 
experimentally.  From  these  analytical  methods,  design  curves  have  been 
produced  for  commonly  used  pipe  sizes  and  fluids.  Similarly,  analytical 
methods  have  been  developed  for  honeycomb  filters  and  calculations  have  been 
made  to  show  shielding  effectiveness  of  this  filter  versus  frequency.  Results 
show  that: 

1.  For  fluids  with  relatively  high  dielectric  constants  and  low  loss 
tangents,  currently  used  designs  can  compromise  the  effectiveness  of  the 
shield. 


2.  A  honeycomb  grating  can  have  its  shielding  effectiveness 
substantially  reduced  if  filled  with  distilled  water  instead  of  air. 

A  problem  encountered  in  this  study  is  that  published  data  on  the 
relative  dielectric  constant  and  loss  tangent  are  incomplete  for  most  fluids 
of  interest. 
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APPENDIX  A: 


WAVE  PROPAGATION  IN  CIRCULAR  METALLIC  WAVEGUIDES 


The  anaLysis  of  circuLar  metallic  waveguides  can  be  found  in  many 
standard  electromagnetic  textbooks.  Significant  results  are  summarized  here. 

A  circuLar  metallic  pipe  supports  both  TE  and  TM  waveguide  modes.  The 
propagation  constants  of  these  modes  are  given  by: 


6  =  /T  k 

n  ,m  r  o 


/  /  c(n,m)  v  2 

\  E  ) 


[Eq  All 


where  =  relative  dielectric  constant  of  the  medium  filling  the  waveguide 

kQ  =  free  space  wave  number 

f  \  =  cutoff  frequency  of  the  TE„  _  or  TM„  „  mode 

cvn,nu  ^  1  n,m  n,m 

f  =  frequency  of  operation. 

Equation  Al  is  valid  for  f  >  fc(n  •  For  TE  modes ,  the  cutoff  frequency  is 
given  by:  ’ 


f 


=  1_ 

c(n,m)  2  it 


n,m  c 


[Eq  A2 1 


where  a  =  radius  of  the  waveguide 

c  =  speed  of  light  in  free  space 

o'  -  the  mth  root  of  Jn'»  the  first  derivative  of  the  nth 
*  order  Bessel  function 

Similarly,  for  TM  modes: 


c(n,m)  2ir  a  /e- 

r 


[Eq  A3] 


where  p  =  the  mth  root  of  J  .  the  nth  order  Bessel  function. 
n,m  n’ 

The  p  and  p'n  mvalues  can  be  found  in  standard  mathematical  tables.’’ 

Tables  Al  aftS  A2  list' values  of  p'  for  m  =  1,  2,  3  and  n  =  0,  1,  2,  3. 

n,m 


4S.  Ramo,  J.  Whinnery,  and  T.  Van  Duzer,  Field  and  Waves  in  Communication 
Electronics  (John  Wiley  and  Sons,  Inc.,  1965);  R.  E.  Collin,  Field  Theory  of 
Guided  Waves  (McGraw-Hill,  Co.,  I960). 

5M.  Abramovitz  and  I.  Stegun,  Handbook  of  Mathematical  Functions  (Dover  Press, 
1970). 
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TabLe  A1 


Roots  of 

J',  Derivative 
n 

of  the 

nth  Order 

Bessel  Function 

n 

m 

0 

1 

2 

3 

i 

3.832 

1.841 

3.054 

4.201 

2 

7.016 

5.331 

6.706 

8.015 

3 

10.173 

8.536 

9.969 

11.346 

Table  A2 

Roots  of  nth  Order 

n 

Bessel  Function 

n 

m 

0 

1 

2 

3 

1 

2.405 

3.832 

5.136 

6.380 

2 

5.520 

7.016 

8.417 

9.761 

3 

8.654 

10.173 

11.620 

13.015 

For  f  <  fc(n  the  (n,m)  waveguide  mode  is  "cutoff."  The  attenuation 
constant  is  given’by: 


a 

n,m 


/F  k 
r  o 


[Eq  A4] 


where  e^,  kQ,  fc(n  and  f  are  as  defined  in  Equation  Al. 

These  results,  typically  derived  for  the  case  in  which  the  waveguide  is 
filled  with  a  lossless  medium,  are  readily  extended  to  the  lossy  case  by 
allowing  to  be  complex.  Denoting  the  complex  dielectric  constant  by  e^*, 

e  *  =  e  -  i—  [Eq  A5] 

r  r  we 

o 

where  e  *  real  relative  dielectric  constant 
r 

o  =  conductivity  of  the  medium 


eQ  =  free~space  permittivity 
w  =  2rf. 
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A  more  convenient  version  of  Equation  A5  is: 


e  *  =  e  (1  -  jtanfi) 
r  r 


[Eq  A6 ] 


where  tan  5  =  — - . 

CO  £  £ 

o  r 

Tan  6  is  known  as  the  loss  tangent  or  dissipation  factor.  The  Loss  tangent  is 
a  measure  of  a  material's  Loss  properties. 

Substituting  e  *  for  £  in  Equation  A1  yieLds  an  expression  for  the 
complex  propagation  constant  S-ja.  The  branches  of  the  complex  square  root 
are  chosen  such  that  a  and  S  are  positive  real.  That  is,  the  soLution 
corresponds  to  a  forward  propagating  wave  with  attenuation. 

It  is  important  to  note  that  there  is  no  real  cutoff  frequency  in  a  Lossy 
waveguide.  Equation  A1  has  a  compLex  soLution  for  all  frequencies,  rather 
than  the  pureLy  reaL  or  purely  imaginary  solutions  present  in  the  Lossless 
case.  The  cutoff  frequency  in  the  lossy  case  can  still  be  defined  as  that  for 
the  lossless  case.  This  definition  is  useful  because  as  the  frequency  changes 
from  f  >  f  to  f  <  f  ,  the  wave  shifts  from  one  dominated  by 

clossless  clossless 

the  propagation  term,  8,  to  one  dominated  by  the  attenuation  term,  a. 
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APPENDIX  B: 


COUPLING  OF  ELECTROMAGNETIC  WAVES  TO  A  CIRCULAR  WAVEGUIDE 


When  an  electromagnetic  wave  is  incident  on  the  open  end  of  a  circular 
waveguide,  some  of  the  wave's  energy  will  be  coupled  into  the  waveguide. 
Analysis  of  the  coupling  factor  is  complicated  by  diffraction  of  the  incident 
wave  by  the  waveguide  end.  Pearson  has  determined  theoretically  the  fields 
within  the  waveguide  for  a  plane  harmonic  electromagnetic  wave  incident  on  a 
semi-infinite  open-ended  perfectly  conducting  circular  waveguide.  This 
problem  has  been  addressed  experimentally  by  Southworth  and  King. 

Experiments  were  done  to  determine  the  directive  properties  of  metal  pipes  and 
horns  when  used  as  electromagnetic  wave  receivers.  This  appendix  presents  the 
basic  textbook,  approach  to  power  transfer  between  transmitter  and  receiver. 

Two  important  antenna  parameters  are  directivity  and  effective 
aperture.  The  directivity  of  an  antenna  is  a  measure  of  the  relative 
concentration  of  radiated  or  received  power  in  different  directions.  It  is 
dependent  on  antenna  type  and  size. 

Effective  aperture  is  defined  as  the  ratio  of  the  power  in  the  antenna 
terminating  impedance  to  the  power  density  of  the  incident  wave.  Under 
conditions  of  maximum  power  transfer  and  zero  antenna  losses,  the  effective 
anerture  is  a  maximum.  Maximum  effective  aperture  is  related  to  directivity 
through  the  expression: 


em 


[Eq  B1  1 


where  D  is  the  directivity,  A„_  is  the  maximum  effective  aperature  and  X  is 

'  em 

the  wavelength. 

Consider  a  transmitting  antenna  with  directivity  D^  and  a  receiving 
antenna  with  maximum  effective  aperture  AgR  (Figure  Bl).  The  antennas  are 
separated  by  a  distance  d.  The  received  power,  PR,  is  given  by: 


P 


R 


P  D 
T  T 


eR 


(Eq  B2  ] 


• 

^J.  D.  Pearson,  "The  Diffraction  of  Electro-Magnetic  Waves  by  a  Semi-Infinite 
Circular  Wave  Guide,"  Proceedings  of  the  Cambridge  Philosophical  Society, 

Vol  49,  Part  4  (1953),  pp  659-667. 

^G.  C.  Southworth,  and  A.  P.  King,  "Metal  Horns  as  Directive  Receivers  of 
Ultra-Short  Waves,"  Proceedings  of  the  I.R.E.,  Vol  27  (February  1939),  • 

pp  95-102. 

°J.  D.  Kraus,  Antennas  (McGraw-Hill,  Co.,  1950);  R.  E.  Collin  and  F.  J. 

Zucker,  Antenna  Theory,  Part  1  (McGraw-Hill,  1969). 
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Figure  Bl.  Transmitting  and  receiving  antennas. 


where  PT  is  the  transmitted  power.  This  formula  assumes  a  LossLess  medium  and 
no  polarization  mismatch.  Using  Equation  Bl,  Equation  B2  can  be  expressed  in 
terms  of  the  directivity  of  the  receiving  antenna,  D^: 


°T  °R  [4 ltd] 


[Eq  B3] 


Therefore,  the  power  transfer  ratio  between  two  antennas  is: 


—  =  D  D  2 

PT  T  R  |_4*d 


[Eq  B4] 


Equation  B4  is  a  far-field  relation  and,  therefore,  will  not  be  valid  if  d  is 
small  compared  to  the  size  of  the  antenna.  The  error  is  small  if: 


d>-H- 


where  1  is  the  maximum  linear  dimension  of  either  antenna. 


[Eq  B5 ] 


The  directivity  of  a  horn  antenna  with  a  large  aperture  can  be  written  in 
terms  of  the  maximum  effective  aperture.  Thus,  from  Equation  Bl: 


D  ■  V  ahx 


[Eq  B6 ] 


where  A_  =  aperture  length  in  tree-space  wavelengths  in  K  plane 
EX 

A  =  aperture  length  in  free-space  wavelengths  in  H  plane 
HX 

Y  =  ratio  of  maximum  effective  aperture  to  physical  aperture. 

An  optimum  horn  has  a  value  of  y  a  0.6.  Thus,  Equation  B6  becomes: 

0  .  7.5  AE1  Ahi  I"  B’l 

The  horns  used  in  the  tests  conducted  for  this  report  have  an  aperture 
size  of  5.85  in.  (14.86  cm)  (E-plane)  by  7.92  in.  (20.12  cm)  (H-plane).  So, 


.  14.86  cm 

aea  =  — a - 

,  20.12  cm 

aha  =  — x - 


[Eq  B8] 


.  • 


Therefore,  the  directivity  of  the  horn  antennas,  Dy,  is  given  by 


Dh  .  7.5 


14.86  20.12 


2242 


[Eq  B9] 


where  X  is  measured  in  cm. 

The  open-ended  waveguide  can  be  considered  as  an  antenna.  The 
directivity  of  an  open-ended  circular  waveguide  is  given  by 


D  3  10.5 
£  X* 


[Eq  B10] 


where  r  is  the  radius  of  the  waveguide.  The  pipes  (circular  waveguides)  used 
in  the  tests  have  radii  of  0.738  in.  (1.874  cm)  and  2.0  in.  (5.08  cm).  Thus, 


D  s  (1.874)2  =  ~  [Eq  811] 

S  X  X 

for  the  small  pipe  and 

0p  *  <5-08>2  =  ^  [Eq  b121 

L  A  A 

for  the  large  pipe.  Again,  A  is  measured  in  cm. 

Figures  B2  and  B3  show  the  placement  of  the  horn  antennas  used  in  these 
tests.  Figure  B2  shows  the  positioning  in  free  space  for  the  reference 
measurements.  The  distance  d  was  measured  from  mouth  to  mouth.  For  testing 
the  waveguide  below  cutoff  theory,  the  configuration  in  Figure  B3  was  used. 
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SHIELDED  ROOM 


Figure  B2.  Experimental  setup  for  reference  measurements. 


SHIELDED  ENCLOSURE 


Figure  B3.  Experimental  setup  for  waveguide  measurements. 


The  waveguide  was  6  in.  (15.24  cm)  tong  with  a  radius  r,  and  d’  is  the 
distance  from  the  horn  mouth  to  the  end  of  the  waveguide. 

As  an  example  calculation  for  coupling  loss,  consider  a.  plane  wave  at  2.5 
GHz  propagating  from  the  transmitting  antenna  toward  the  large  pipe.  Some  of 
the  transmitted  wave  power  will  be  coupled  into  the  pipe  and  some  will  be 
diffracted,  which  results  in  a  coupling  loss.  Some  of  the  wave  power  that 
coupled  into  the  pipe  will  retransmit  at  the  other  end  of  the  pipe.  Again, 
there  will  be  a  coupling  loss.  Finally,  the  resultant  wave  power  will  be 
sampled  by  the  receiving  horn. 

The  coupling  loss  is  one-half  the  dB  difference  between  horn-to-horn 
power  transfer  and  horn-to-pipe-to-horn  power  transfer.  If  ^h-H  den°t-es  the 
power  transfer  ratio  between  the  two  horns  in  free  space  and  P^_p_^  denotes 
the  power  transfer  ratio  between  the  two  horns  separated  by  the  pipe,  then 
coupling  loss,  Lc,  s: 


P  -  P 

H-H  H-P-H  jnX 

L  =  - 2 -  (*n  dB' 


[Eq  B13] 


Let  the  horn-to-horn  distance  d  (Figure  B2)  be  78  in.  (198  cm).  Then: 

P„ 


p  =  —  = 

H-H  PT 


°T  DR  j^TTdj 

*  °K  D«  [ra]  2 


Using  Equation  B9: 


Recalling: 


Thus , 


H-H  L 


Expressed  in  dB: 


P  s 

2242 

2 

X 

H-H 

_x2  J 

4xd 

10 

X  =  ■§■  - 

3  x  10 

=  12 

f  2 

.5  x 

109 

C  22421  2 

r 

12 

1 

|_!44  J 

1  4tt 

( 198)  J 

;-H  *  10  lo8l0 

(.00564) 

=  .00564 


[Eq  B14] 


[Eq  B15] 


[Eq  B16] 


[Eq  B17 ] 


[Eq  B18 ] 


35 


Let  the  horn-to-pipe  distance  d'  (Figure  B2b)  be  36  in.  (91.4  cm).  Then: 


ph-?l  “  °H  °P  [4^] 

Using  Equations  B9  and  B12: 


P  _f~2242~l  [~851 

H-P.  ,2  2 

L  L  X  J  LX  J 


Thus, 


p  _  |~2242 
H-Pl  |_144 


Expressed  in  dB: 

ph-pl  =  -20  dB 


85 1] 

144]  [j 


12 


4n  (91.4) 


=  .01 


Because  of  reciprocity: 


PH-pL  =  PpL"H 


So: 


ph-pl-h  =  2  ph-pl  =  "40  dB 

Therefore, 


L 


C 


-22.5  -  (-40) 
2 


=  8.75  dB 


[Eq  B19 1 


[Eq  B20 ] 


[Eq  B21] 


[Eq  B22 ] 


[Eq  B23 ] 


[Eq  B24J 


[Eq  B25 ] 


With  the  transmitter  supplying  5  W  at  2.5  GHz,  a  value  of  +5  dBm  was 
measured  by  the  receiver  for  the  reference  measurement.  The  waveguide 
measurement  was  -10  dBm.  This  resuLts  in  a  measured  coupling  loss  of: 

Lc  =  =  7.5  dB  [Eq  B26 ] 

The  error  in  reading  the  receiver  is  +1/2  dB,  giving  a  +3/4  dB  error  in 
the  measured  coupling  loss.  Greater  error  probably  exists  in  the  reference 
measurement.  Because  of  cable  constraints  (too  short),  it  was  not  possible  to 
move  the  receiving  horn  farther  than  36  in.  (91.4  cm)  from  the  shielded 
enclosure.  This  put  the  receiving  horn  in  the  "standing  wave"  field  between 
the  enclosure  and  the  transmitting  horn,  resulting  in  a  non-f ree-f ield 
measurement  for  the  reference  reading. 


APPENDIX  C: 


FORTRAN  PROGRAMS  TO  CALCULATE  WAVEGUIDE  ATTENUATION 
RESULTING  FROM  DIELECTRIC  LOADING 


Programs  Pipe  1A  and  Pipe  IB  (Figures  Cl  and  C2)  calculate  the  overall 
attenuation  in  dB  of  an  electromagnetic  wave  propagating  through  a  length  of 
circular  waveguide  (pipe).  Program  Pipe  lA  uses  one  relative  dielectric 
constant  and  loss  tangent  to  compute  the  attenuation  for  a  band  of 
frequencies.  Program  Pipe  IB  uses  a  separate  relative  dielectric  constant  and 
loss  tangent  for  each  frequency. 

The  attenuation  factor  (in  Np/cm)  of  a  circular  waveguide  is  the  real 
part  of  the  complex  propagation  constant. 


Y 


2 

u  pe 


[Eq  Cl] 


where  u>  =  angular  frequency  (rad/sec) 

u  =  permeability  of  the  material  in  the  waveguide 
e  =  permittivity  of  the  material  in  the  waveguide,  and: 

(ha) 

h  «  - [Eq  C2 ] 

nm  a 

where  a  is  the  pipe  radius  and  (ha)nm  is  the  mth  root  of  a  Bessel  function  of 

the  nth  order.  Values  of  (ha)  for  TM  modes  are  given  in  Appendix  A,  Table 

A2;  values  of  (ha)„_  for  TE  mooes  are  given  in  Table  Al. 
nm 

The  cutoff  frequency  of  the  (n,m)  mode  is  the  frequency  for  which  the 
real  part  of  y  is  zero.  Attenuation  through  the  waveguide  is  a  minimum  at 
the  cutoff  frequency. 


Since 


[Eq  C3 ] 


c  = 


/rr 

O  Q 


[Eq  C4j 


where  c  is  the  speed  of  light  in  a  vacuum. 
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PROGRAM  PIPE1 A ( INPUT , TAPE** INPUT ) 

A  PROGRAM  TO  CALCULATE  THE  ATTENUATION  IN  A  CIRCULAR  METALLIC 
UAVE6UIDE  FILLED  UITH  LOSSY  DIELECTRIC. 

FEATURES:  WAVEGUIDE  ATTENUATION  CALCULATIONS 
CALCULATIONS  OVER  A  FREQUENCY  RAN6E 


ELECTROMAGNETIC  COUPLING  THROUGH  A  FLUID  FILLED  PIPE  IS  MODELLED 
BY  THE  ATTENUATION  IN  THE  TEH1)  CIRCULAR  UAVE6UIDE  MODE.  THIS 
REPRESENTS  THE  FOLLOWING  TWO  CASES. 

1)  WORST-CASE  ANALYSIS  -  ALL  THE  ELECTROMAGNETIC  ENERGY 
INCIDENT  ON  THE  C0UPLIN6  RE6I0N  OF  THE  PIPE  IS  COUPLED  INTO 
THE  LOWEST  ORDER  WAVEGUIDE  MODE. 

2)  ACCURATE  KNOWLEDGE  OF  THE  EXCITATION  MECHANISM  WHEREBY 
ELECTROMAGNETIC  ENER6Y  IS  COUPLED  INTO  THE  WAVEGUIDE  IS 
UNAVAILABLE,  MAKING  A  MORE  ACCURATE  MODAL  ANALYSIS 
IMPOSSIBLE. 


USAGE:  COMPILE  UITH  THE  FTN5  (FORTRAN  77)  COMPILER. 

THE  PROGRAM  IS  WRITTEN  TO  BE  USED  INTERACTIVELY.  INPUT 
IS  ENTERED  FROM  THE  TERMINAL  IN  RESPONSE  TO  QUESTIONS. 
OUTPUT  IS  PRINTED  TO  THE  TERMINAL. 

FREQUENCY  AND  OVERALL  ATTENUATION  IS  ALSO  WRITTEN  TO  AN 
OUTPUT  FILE  USING  THE  FORMAT.  THIS  FILE  CAN  THEN  BE 
PRINTED  OR  PLOTTED  AT  A  LATER  TIME.  THE  USER  IS 
PROMPTED  FOR  AN  OUTPUT  FILE  NAME.  A  DEFAULT  NAME 
COUTPF")  IS  PROVIDED  IF  THE  USER  ENTERS  A  CARRIAGE 
RETURN  IN  RESPONSE  TO  THIS  QUESTION.  THE  OUTPUT  FILE  IS 
AUTOMATICALLY  SAVED  UNLESS  THE  DEFAULT  NAME  IS  USED. 

WHEN  USIN6  THE  DEFAULT  NAME,  THE  OUTPUT  FILE  REMAINS 
LOCAL  AND  MUST  BE  EXPLICITLY  SAVED  (IF  DESIRED)  USIN6 
SYSTEM  LEVEL  COMMANDS. 

INPUT:  PIPE  RADIUS  IN  CM 
PIPE  LEN6TH  IN  CM 
DIELECTRIC  CONSTANT 
LOSS  TANGENT 

MINIMUM  FREQUENCY  IN  GHZ 


Figure  Cl.  Program  Pipe  1A 
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I 


46  C 

47  C 

48  C 

49  C 

50  C 
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54  C 
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59  C 

60 
61 
62 

63 

64 

65 
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67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87  10 

88 

89 

90 


MAXIMUM  FREQUENCY  IN  GHZ 
FREQUENCY  INCREMENT  IN  GHZ 
OUTPUT  FILE  NAME 

OUTPUT:  ATTENUATION  CONSTANT  (ALPHA)  IN  NEPERS/CM 
ATTENUATION  CONSTANT  (ALPHA)  IN  DB/CH 
OVERALL  ATTENUATION  IN  DB 

OUTPUT  TO  FILE:  FREQUENCY  IN  GHZ  VS.  OVERALL  ATTENUATION  IN  OB 

LANGUAGE:  FORTRAN  77 

PROGRAMMER:  SCOTT  L.  RAY 
DATE:  7/11/83 


COMPLEX  GAMMA, J 
REAL  LENGTH,KO,KC 
CHARACTER  *7  NANE.NTEMP 
PI»3. 14159265 
J»(0. 0,1.0) 

P1 1 P=1 - 84118 
A1=AL0G10(2. 71828183) 

PRINT  1000 
READ  ♦, RADIUS 
KC*P1 IP/RADIUS 
PRINT  1010 
READ  *, LENGTH 
PRINT  1020 
READ  *  EPSR 

FC=1.0/(2.0*PI)*(29.997925/SQRT(EPSR))*KC 

PRINT  1030 

READ  «,DISP 

PRINT  1040 

READ  *,FMIN 

PRINT  1050 

READ  *,FHAX 

PRINT  1060 

READ  *,FINC 

NAMEs'OUTPF' 

PRINT  1070 

READ ( 6 , 2030 , END= 1 0 )  NTEMP 

NAME*NTEMP 

REWIND  6 

OPEN (3, FILE “NAME) 

PRINT  2060 
PRINT  2040 


Figure  Cl.  (Cont'd). 


39 


91  PRINT  1080,  rtDIUS, LENGTH 

92  PRINT  1090,EPSR,DISP 

93  PRINT  2000, FC 

94  PRINT  2010 

95  PRINT  2050 

94  DO  20  FR£Q*FMIN, (FHAX+FINC/2.0) ,FINC 

97  K0s2.0*PI*FREQ/30.0 

98  GAHMA=CSQRT(KC*KC-EPSR*(1 .0-J*DISP)*KO*KO) 

99  ALPHA*REAL(GAMMA) 

100  ALPDBS20.0*A1* ALPHA 

101  ATTEN=ALPDB*LEN6TH 

102  PRINT  2020, FREQ, ALPHA, ALPDB,ATTEN 

103  URITE<3,*)  FREQ, ATTEN 

104  20  CONTINUE 

105  PRINT  2040 

104  CLOSE (3) 

107  IF(NAME.NE.'OUTPF')  CALL  PF(  'SAVE' , NAME) 

108  STOP 

109  1000  FORMAT ('ENTER  PIPE  RADIUS  IN  CM') 

110  1010  FORMAT ('ENTER  PIPE  LENGTH  IN  CM') 

111  1020  FORMAT ( 'ENTER  DIELECTRIC  CONSTANT') 

112  1030  FORMAT ('ENTER  THE  LOSS  TANGENT') 

113  1040  FORMAT ('ENTER  MINIMUM  FREQUENCY  IN  GHZ') 

114  1050  FORMAT ('ENTER  MAXIMUM  FREQUENCY  IN  GHZ') 

115  1040  FORMAT ('ENTER  FREQUENCY  INCREMENT  IN  GHZ') 

114  1070  FORMAT ('ENTER  THE  OUPUT  FILE  NAME,  (DEFAULT  OUTPF)') 

117  1080  FORMAT (/'RADIUS  =  ',F7.3,'  CM' ,5X, 'LENGTH  *  ',F8.3,'  CM') 

118  1090  FORMAT( 'EPSILON  *  ' ,F4.2,8X, 'LOSS  TANGENT  =  ',F4.4> 


119  2000  FORMAT < / '  CUTOFF  FREQUENCY  =  ',F7.3,'  GHZ') 

120  2010  F0RMAT(/3X,'FREQ  (GHZ)', 3X, 'ALPHA  (NEPERS/CM)' ,3X,' ALPHA  (DB/CM) 

121  *',3X, 'ATTEN  (DB)') 

122  2020  F0RHAT(5X,F5.2,9X,1PE10.3,10X,0PF5.2,8X,F7.2> 

123  2030  F0RMAT(A7) 

124  2040  FORMATt' . . 

1 25  * - ') 

124  2050  F0RMAT(3X,' . - . . 

127  *- . ') 

128  2040  FORMAT ( // ) 

129  END 


Figure  Cl.  (Cont'd). 
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39  C 

40  C 
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PROGRAM  PIPE 1B< INPUT, TAPE6*I NPUT ) 

A  PROGRAM  TO  CALCULATE  THE  ATTENUATION  IN  A  CIRCULAR  METALLIC 
UAVE6UIDE  FILLED  WITH  LOSSY  DIELECTRIC. 

FEATURES:  UAVE6UIDE  ATTENUATION  CALCULATIONS 
USER  SPECIFIES  EACH  FRE8UENCY 
ABILITY  TO  CHANGE  MATERIAL  PARAMETERS  AT  EACH  NEU 
FREQUENCY 


ELECTROMAGNETIC  COUPLING  THROUGH  A  FLUID  FILLED  PIPE  IS  MODELLED 
BY  THE  ATTENUATION  IN  THE  TE( 1 1 )  CIRCULAR  UAVE6U1DE  MODE.  THIS 
REPRESENTS  THE  FOLLOWING  TWO  CASES. 

1)  WORST-CASE  ANALYSIS  -  ALL  THE  ELECTROMAGNETIC  ENERGY 
INCIDENT  ON  THE  COUPLING  REGION  OF  THE  PIPE  IS  COUPLED  INTO 
THE  LOWEST  ORDER  WAVEGUIDE  MODE. 

2)  ACCURATE  KNOWLEDGE  OF  THE  EXCITATION  MECHANISM  WHEREBY 
ELECTROMAGNETIC  ENER6Y  IS  COUPLED  INTO  THE  WAVEGUIDE  IS 
UNAVAILABLE,  MAKING  A  MORE  ACCURATE  NODAL  ANALYSIS 
IMPOSSIBLE. 


USAGE:  COMPILE  WITH  THE  FTNS  (FORTRAN  77)  COMPILER. 

THE  PROGRAM  IS  WRITTEN  TO  BE  USED  INTERACTIVELY.  INPUT 
IS  ENTERED  FROM  THE  TERMINAL  IN  RESPONSE  TO  QUESTIONS. 
OUTPUT  IS  PRINTED  TO  THE  TERMINAL. 

WHEN  CALCULATIONS  ARE  COMPLETED  FOR  ONE  FREQUENCY  AND 
SET  OF  MATERIAL  PARAMETERS,  THE  USER  IS  PROMPTED  FOR 
NEXT  FREQUENCY,  DIELECTRIC  CONSTANT,  AND  LOSS  TANGENT. 
PREVIOUS  VALUES  ARE  USED  AS  DEFAULTS.  A  DEFAULT  VALUE 
IS  OBTAINED  BY  HITTING  CARRIAGE  RETURN  IN  RESPONSE  TO 
THE  PROMPTING  QUESTION. 

FREQUENCY  AND  OVERALL  ATTENUATION  ARE  WRITTEN  TO  AN 
OUTPUT  FILE  USING  THE  FORMAT.  THIS  FILE  CAN  THEN  BE 
PRINTED  OR  PLOTTED  AT  A  LATER  TIME.  THE  USER  IS 
PROMPTED  FOR  AN  OUTPUT  FILE  NAME.  A  DEFAULT  NAME 
("OUTPF")  IS  PROVIDED  IF  THE  USER  ENTERS  A  CARRIAGE 
RETURN  IN  RESPONSE  TO  THIS  QUESTION.  THE  OUTPUT  FILE  IS 
AUTOMATICALLY  SAVED  UNLESS  THE  DEFAULT  NAME  IS  USED. 


Figure  C2.  Program  Pipe  IB. 
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46  C 

47  C 

48  C 
4?  C 

50  C 

51  C 

52  C 

53  C 

54  C 

55  C 

56  C 

57  C 

58  C 

59  C 

60  C 

61  C 

62  C 

63  C 

64  C 

65  C 

66  C 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85  10 

86 

87 

88 

89 

90  20 


UHEN  USING  THE  DEFAULT  NAME,  THE  OUTPUT  FILE  REMAINS 
LOCAL  AND  MUST  BE  EXPLICITLY  SAVED  (IF  DESIRED)  USING 
SYSTEH  LEVEL  COMMANDS. 

INPUT:  PIPE  RADIUS  IN  CM 

PIPE  LENGTH  IN  CM 

DIELECTRIC  CONSTANT 
LOSS  TANGENT 
FREOUENCY  IN  GHZ 
OUTPUT  FILE  NAME 

OUTPUT:  ATTENUATION  CONSTANT  (ALPHA)  IN  NEPERS/CM 

ATTENUATION  CONSTANT  (ALPHA)  IN  DB/CM 

OVERALL  ATTENUATION  IN  DB 

OUTPUT  TO  FILE:  FREQUENCY  IN  GHZ  VS.  OVERALL  ATTENUATION  IN  DB 

LAN6UA6E:  FORTRAN  77 

PROGRAMMER:  SCOTT  L.  RAY 
DATE:  7/11/83 


REAL  LENGTH, KQ,KC 
COMPLEX  J, GAMMA 
CHARACTER  *7  NAME ,NTEMP 
PI=3. 14159265 
J*(0. 0,1.0) 

P11P*1 .84118 

A1=AL0610(2. 71828183) 

EPSR*1 .0 

DISP=0.0 

PRINT  1000 

READ  *, RADIUS 

KC=P1 IP/RADIUS 

PRINT  1010 

READ  *, LENGTH 

NAHE='OUTPF' 

PRINT  1020 

READ(6,2010,ENDs10)  NTEMP 

NAME*NTEMP 

REUIND  6 

OPEN  (3,  FILENAME) 

PRINT  2020 
PRINT  1030 
READ  «,FREQ 
REUIND  6 


Figure  C2.  (Cont ’d) 
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91 

KO=2.0*PI*FREQ/30.0 

92 

PRINT  1040,  EPSR 

93 

RE AD ( 6 , * , END=30 ) TEHP 

94 

EPSR=TEMP 

95 

30 

REUIND  6 

96 

FC=1 .0/(2. 0*PI ) *  <  29 . 997925/SQRT (EPSR) )*KC 

97 

PRINT  1050, DISP 

98 

READ(6,*,END«40)  TEMP 

99 

DISPaTEMP 

100 

40 

REUIND  6 

101 

GAMMAaCSQRKKC*KC-EPSR*< 1 .0-J*DISP)*K0*K0> 

102 

ALPHAaRE AL ( 6ANNA ) 

103 

ALPDI*20.0*A1*ALPHA 

104 

ATTEN=ALPDB*LENGTH 

105 

PRINT  1060, ALPHA, ALPDB 

106 

PRINT  1065,ATTEN 

107 

URITE < 3 , ♦ )  FREQ,ATTEN 

108 

PRINT  2020 

109 

TEMPaFREQ 

110 

PRINT  1070, FREQ 

111 

PRINT  1080 

112 

READ(6,*,ENDa20)  TEMP 

113 

FREQaTEMP 

114 

IF(FREQ.NE.O.O)  GO  TO  20 

1  i5 

CLOSE  <  3 ) 

116 

IF < NAME .NE. "OUTPF ' )  CALL  PF( 'SAVE' , NAME) 

117 

STOP 

118 

1000 

FORNAT( 'ENTER  PIPE  RADIUS  IN  CH') 

119 

1010 

FORMAT ('ENTER  PIPE  LENGTH  IN  CM') 

120 

1020 

FORMAT ('ENTER  THE  OUPUT  FILE  NAME,  (DEFAULT  OUTPF)') 

121 

1030 

FORMAK 'ENTER  THE  FREQUENCY  IN  GHZ' ) 

122 

1040 

FORMAT ('ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  ',F6.2) 

123 

1050 

FORMAK 'ENTER  THE  LOSS  TANGENT,  DEFAULT  ',F5.3) 

124 

1060 

FORMAK//'  ALPHA  a  ',1P1E10.3,'  NEPERS/CM  <',0PF5.2,'  DB)') 

125 

1065 

FORMAT ('  TOTAL  ATTENUATION  =  ',F7.2,'  DB') 

126 

1070 

FORMAK 'ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  ',F8.3) 

127 

1080 

FORMAK 'ENTER  0  TO  QUIT') 

128 

2010 

FORMAT (A7) 

129 

2020 

FORMAK//' . ') 

130 

END 

k» 


I 


Figure  C2.  (Cont'd)* 
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[Eq  C5] 


f  =  [Eq  C5] 

c  2  w/e 

r 

where  e  is  the  relative  dielectric  constant^  of  the  material  in  the 
waveguiSe.  Because  only  the  real  part  of  y  is  used  to  define  cutoff,  the 
loss  tangent  (the  imaginary  part  of  the  complex  permittivity)  is  neglected. 

Y  is  computed  at  line  number  98  in  program  Pipe  1A  (line  number  101  in 
program  Pipe  IB).  The  program  variables  used  in  the  computation  of  y  are 
defined  in  Table  Cl.  Both  programs  use  the  same  variables,  so  only  one  table 
is  presented.  The  overall  attenuation  in  dB  resulting  from  the  length  of  pipe 
is  computed  at  line  number  101  (104)  in  program  Pipe  1A  (Pipe  IB). 

A  sample  run  of  each  program  for  distilled  water  is  presented  in  Figures 
C3  and  C4.  A  pipe  of  radius  1.0  in.  (2.54  cm)  and  length  6  in.  (15.24  cm)  was 
modeled.  Table  C2  lists  the  relative  dielectric  constants  and  loss  tangents 
used  for  the  frequency  band  of  1.0  to  10.0  GHz.  These  values  were 
interpolated  linearly  from  data  given  by  von  Hippel  for  300  MHz,  3  GHz,  and  10 
GHz. ^  Because  Program  Pipe  1A  uses  one  relative  dielectric  constant  and  one 
loss  tangent  for  a  band  of  frequencies,  the  average  value  of  these  parameters 
over  the  frequency  band  of  1.0  GHz  to  10  GHz  was  used. 

Table  C3  is  a  summary  of  the  attenuation  output  data  generated  by  these 
programs.  Note  that  between  1  and  5  GHz,  Pipe  IB  predicts  less  attenuation 
than  Pipe  1A  and  vice-versa  between  6  and  10  GHz.  This  occurs  because 
attenuation  is  very  dependent  on  loss  tangent.  Between  1  and  5  GHz,  the  loss 
tangent  used  by  Pipe  1A  is  greater  than  that  used  by  Pipe  IB,  resulting  in 
greater  attenuation.  Between  6  and  10  GHz,  the  loss  tangent  used  by  Pipe  1A 
is  less  than  that  used  by  Pipe  IB,  resulting  in  less  attenuation. 


Q 

■'A.  R.  von  Hippel. 
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Table  Cl 


Variables  Used  in  Programs  Pipe  1A  and  Pipe  IB 


/ariable 

Definition 

Line  nos. 
where  defined 
(Pipe  IA,  Pipe  IB) 

A1 

L°8l0  e 

66,73 

DISP 

tan  4,  Loss  tangenc 

76,98 

EPSR 

e  ,  Relative  dielectric  constant 
r 

73,93 

FREQ 

f,  Frequency  in  CHz 

96,(89  or  112) 

J 

j. 

64,71 

LENGTH 

1,  Length  of  pipe  in  cm 

71,80 

PI 

* 

63,70 

P11P 

(ha)nnl  for  TE^  mode 

65,72 

RADIUS 

r,  Radius  of  pipe  in  cm 

68,77 

KC 

kc»  (ha)^/r  (r  in  cm) 

69,78 

KO 

k„,  2»f/30  (f  in  CHs) 

97,91 

GAMMA 

Tn/k  *  -  t  (1  -  jtani)  k  * 
c  r  0 

98,101 

ALPHA 

a,  Real  part  of  y 

99,102 

ALPDB 

dB  a 

Attenuation  in - ,20  Log  1Ae 

CfB 9  1U 

100,103 

ATTEN 

Overall  attenuation  in  dB,  ALPDB  x  1 

101,104 

FC 

,  1  29.997925, 

Cutoff  frequency,  ^  — /f= - k c 

r 

74,96 

Table  C2 

Relative  Dielectric  Constants  and  Loss  Tangents  for 
Distilled  Water  at  25°C 


frequency  Relative  dielectric  constant 
(GHz )  Ur> 


toss  tangent 
(tan  5) 


1.0 

77.3 

.052 

2.0 

77.0 

.105 

3.0 

76.7 

.157 

4.0 

73.6 

.212 

5.0 

70.5 

.266 

6.0 

67.4 

.321 

7.0 

64.3 

.376 

8.0 

61.2 

.430 

9.0 

58.1 

.485 

10.0 

55.0 

.540 

cf  «  68.11 


tan6  ■  .294 
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TERMINAL:  140 

in  00 

U  OF  I  CYBER  175-175  ACCOUNTS  ONLY  NOS  1.4  -  501/498. 

SIGMON:  3KEMVUJ 
PASSWORD 

ilium 

TERMINAL:  140,  TTY 

RECOVER/  CHARGE:  -CERL 

LAST  RECORDED  SIGNON  AT  16:03  08/24/83 
■**  BEGIN  CERL  BILLING  ** 

/GET,PIPE1A 

/R.FTN5,I=PIPE1A,L=0 

0.147  CP  SECONDS  COMPILATION  TIME. 

/LGO 

ENTER  PIPE  RADIUS  IN  CM 
?  2.54 

ENTER  PIPE  LENGTH  IN  CM 
?  15.24 

ENTER  DIELECTRIC  CONSTANT 

?  68.11 

ENTER  THE  LOSS  TANGENT 
?  .294 

ENTER  NIN1HUN  FREQUENCY  IN  GHZ 

?  1.0 

ENTER  MAXIMUM  FREQUENCY  IN  GHZ 

?  10.0 

ENTER  FREQUENCY  INCREMENT  IN  GHZ 

?  1.0 

ENTER  THE  OUPUT  FILE  NAME,  (DEFAULT  OUTPF) 

?  UATER 


RADIUS  =  2.540  CM 

EPSILON  =  68.11 

CUTOFF  FREQUENCY  = 


LENGTH  =  15.240  CM 

LOSS  TANGENT  =  .2940 

.419  GHZ 


FREQ  (GHZ) 

ALPHA  (NEPERS/CM) 

ALPHA  (DB/CM) 

ATTEN  (DB) 

1.00 

2.757E-01 

2.39 

36.49 

2.00 

5.138E-01 

4.46 

68.02 

3.00 

7.615E-01 

6.61 

100.80 

4.00 

1.011E+00 

8.78 

133.84 

5.00 

1 .261 E+00 

10.96 

166.98 

6.00 

1.512E+00 

13.13 

200.17 

7.00 

1 . 763E+00 

15.31 

233.39 

8.00 

2.01 4E+00 

17.49 

266.62 

9.00 

2.265E+00 

19.68 

299.87 

10.00 

2.51 7E+00 

21.86 

333.12 

0.036  CP  SECONDS  EXECUTION  TIME. 


Figure  C3.  Sample  run  of  program  Pipe  1A  for  distilled  water 


TERMINAL:  113 

83/08/25.  14.37.56. 

U  OF  I  CYBER  175-175  ACCOUNTS  ONLY  NOS  1.4  -  501/498. 


SIGNON:  3KEMVUJ  • 

PASSWORD 

MINIM 

TERMINAL:  113,  TTY 

RECOVER/  CHARGE:  -CERL 

LAST  RECORDED  SIGNON  AT  14:36  08/25/83  • 

*•  BEGIN  CERL  BILLING  •* 

/GET ,P1PE1B 
/R.FTN5,I=P1PE1B,L=0 

0.134  CP  SECONDS  COMPILATION  TIME. 

/LGO 

ENTER  PIPE  RADIUS  IN  CM  • 

?  2.54 

ENTER  PIPE  LENGTH  IN  CM 
?  15.24 

ENTER  THE  OUPUT  FILE  NAME,  (DEFAULT  DUTPF) 

?  UATER2 

.  • 


ENTER  THE  FREQUENCY  IN  GHZ 

?  1.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  1.00 
?  77.3 

ENTER  THE  LOSS  TANGENT,  DEFAULT  0.000 
?  .052 


ALPHA  =  5.206E-02  NEPERS/CM  (  .45  DB> 

TOTAL  ATTENUATION  =  6.89  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  1.000 
ENTER  0  TO  OUIT 
?  2.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  77.30 
?  77.0 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .052 
?  .105 


ALPHA  =  1.965E-01  NEPERS/CM  <  1.71  DB) 

TOTAL  ATTENUATION  =  26.02  DB 

Figure  C4.  SampLe  run  of  Program  Pipe  IB  for  distilled  water. 
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ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  2.000 
ENTER  0  TO  QUIT 
?  3.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  77.00 
?  76.7 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .105 
?  .157 


ALPHA  =  4.344E-01  NEPERS/CH  (  3.77  DB) 
TOTAL  ATTENUATION  =  57.50  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  3.000 
ENTER  0  TO  QUIT 
?  4.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  76.70 
?  73.6 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .157 

?  .212 


ALPHA  =  7.614E-01  NEPERS/CK  (  6.61  DB) 
TOTAL  ATTENUATION  -  100.79  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  4.000 
ENTER  0  TO  QUIT 
?  5.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  73.60 
?  70.5 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .212 
?  .266 


ALPHA  =  1 . 163E+00  NEPERS/CH  (10.10  DB) 

TOTAL  ATTENUATION  =  153.98  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  5.000 
ENTER  0  TO  QUIT 
?  6.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  70.50 
?  67.4 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .266 
?  .321 


Figure  C4.  (Cont'd) 


ALPHA  =  1.639E+00  NEPERS/CH  (14.24  DB) 

TOTAL  ATTENUATION  =  216.99  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  6.000 
ENTER  0  TO  QUIT 
?  7.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  67.40 
?  64.3 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .321 
?  .376 


ALPHA  s  2.177E+00  NEPERS/CM  (18.91  DB) 
TOTAL  ATTENUATION  *  288.20  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  7.000 
ENTER  0  TO  QUIT 
?  8.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  64.30 
?  61.2 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .376 
?  .43 


ALPHA  =  2.762E+00  NEPERS/Ch  (23.99  DB) 

TOTAL  ATTENUATION  *  365.57  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  8.000 
ENTER  0  TO  QUIT 
1  9.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  61.20 
?  38.1 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .430 
?  .485 


Figure  C4.  (Cont'd) 


1 
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ALPHA  =  3.395^*00  NEPERS/CH  (29.49  DB) 
TOTAL  ATTENUATION  =  449.39  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  9.000 
ENTER  0  TO  QUIT 
?  10.0 

ENTER  THE  DIELECTRIC  CONSTANT,  DEFAULT  58.10 
?  55.0 

ENTER  THE  LOSS  TANGENT,  DEFAULT  .485 
?  .54 


ALPHA  =  4.061 E +00  NEPERS/CH  (35.28  DB) 
TOTAL  ATTENUATION  =  537.63  DB 


ENTER  THE  FREQUENCY  IN  GHZ,  DEFAULT  10.000 
ENTER  0  TO  QUIT 
?  0 

0.062  CP  SECONDS  EXECUTION  TINE. 

/BYE 

3KEHVUJ  COSTS!  61.219  SRUS  AT  1.0059  =  9.36 


Figure  CA.  (Cont'd). 
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Table  C3 


Summary  of  Output  From  Sample  Runs* 


Atten  (dB) 

Frequency  Pipe  1A  Pipe  IB 

(GHz)  (Constant  e  and  tan  6)  (Variable  e  and  tan  6) 

r  r 


1.0 

36.49 

6.89 

2.0 

68.02 

26.02 

3.0 

100.80 

57.50 

4.0 

133.84 

100.79 

5.0 

166.98 

153.98 

6.0 

200.17 

216.99 

7.0 

233.39 

288.20 

8.0 

266.62 

365.57 

9.0 

299.87 

449.39 

10.0 

333.12 

537.63 

*See  Figures  C3  and  C4 


APPENDIX  D: 


DETERMINATION  OF  MATERIAL  PARAMETERS 


There  are  few  published  measurements  of  the  relative  dielectric 
constant,  e  ,  and  the  loss  tangent,  tan  5,  in  the  microwave  bands  for 
fluids.  Von  Hippel  provides  data  at  discrete,  widely  spaced  frequencies.10 
Since  the  relative  dielectric  constant  varies  slowly  with  frequency  and  does 
not  have  a  great  effect  on  the  overall  waveguide  attenuation,  an  interpolated 
value  can  be  used  for  e  .  However,  the  same  cannot  be  said  for  the  loss 
tangent.  r 

Accurate  values  of  e  ,  and  tan  6  are  required  to  realistically  compare 
theory  with  experiment,  fn  the  microwave  frequency  band,  loss  tangents  can 
vary  widely  with  frequency  and  from  sample  to  sample,  depending  on  the  purity 
of  the  sample.  Most  importantly,  an  accurate  value  for  the  loss  tangent  is 
needed.  As  shown  in  Chapter  3,  a  small  variation  in  tan  6  can  result  in  a 
large  change  in  waveguide  attenuation. 

Loss  tangents  were  calculated  using  transmission  measurements  through  a 
dielectric  slab.  The  dielectric  slab  consisted  of  a  plexiglass  tank,  filled 
with  the  desired  fluid.  The  interior  thickness  of  the  tank  was  6  in.  (15.24 
cm).  The  walls  consisted  of  0.5-in.  ( 1 .27-cm)-thick  plexiglass. 

A  transmitter  and  waveguide  horn  antenna  were  placed  on  one  side  of  the 
tank,  and  a  receiver  and  another  waveguide  horn  antenna  on  the  other  side. 
Precautions  such  as  placing  Che  tank  in  a  shielded  enclosure  wall  were  taken 
to  reduce  multipath  effects. 

Transmission  measurements  were  referenced  to  free-space  measurements  over 
the  same  distance.  This  eliminated  the  error  caused  by  power  variation  with 
frequency  and  minimized  the  error  caused  by  variations  in  the  waveguide  horns' 
efficiency  with  frequency. 

Knowing  the  relative  dielectric  constants  of  the  plexiglass  and  fluid, 
the  loss  tangent  of  the  plexiglass,  and  the  thickness  of  the  plexiglass  and 
fluid,  the  loss  tangent  of  che  fluid  from  a  transmitted  power  measurement  can 
be  calculated.  The  equations  derived  in  Appendix  F  are  used  for  this 
computation.  These  calculations  were  performed  for  three  fluids:  carbon 
tetrachloride,  methanol,  and  distilled  water. 

In  these  calculations,  the  relative  dielectric  constant  for  plexiglass 
was  assumed  to  be  2.5.  The  loss  tangent  for  plexiglass  was  assumed  to  be  0. 
Relative  dielectric  constants  of  the  fluids  were  interpolated  from 
published  measurements.  1 


|®A.  R.  von  Hippel 
*A.  R.  von  Hippel 
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Tables  D1  through  D3  list  the  measured  values  of  tan  6  along  with  the 
published  values.  The  measured  values  of  tan  6  were  consistently  lower  than 
the  published  values,  except  for  carbon  tetrachloride,  for  which  the  loss 
tangent  was  too  small  to  measure  by  the  above  method.  The  discrepancy  is 
probably  explained  by  the  nonideal  conditions  of  the  experiment.  The 
theoretical  analysis  was  based  on  plane  wave  propagation  and  infinite  9labs 
whereas  the  measurements  were  performed  with  horn  antennas  and  a  finite  slab. 

To  estimate  a  more  accurate  value  of  the  loss  tangent,  the  measured 
values  of  tan  6  were  used  to  show  the  trend  (shape  of  the  curve)  of  the 
data.  The  published  data  provide  a  reference  level.  The  last  column  of 
Tables  D1  and  D2  contains  the  corrected  values  of  loss  tangents.  For  carbon 
tetrachloride,  the  published  value  was  used  throughout  the  band. 


Table  D1 


Relative  Dielectric  Constant  and  Loss  Tangent  Data 
for  Carbon  Tetrachloride* 


Frequency 


tan  6 


(GH7.) 

e 

r 

Published*** 

Measured** 

Corrected 

2.5 

2.17 

0.0004 

3.0 

2.17 

0.0004 

0.0004 

3.5 

2.17 

0.0004 

4.0 

2.17 

0.0004 

4.5 

2.17 

0.0004 

5.0 

2.17 

0.0004 

5.5 

2.17 

0.0004 

6.0 

2.17 

0.0004 

6.5 

2.17 

0.0004 

7.0 

2.17 

0.0004 

7.5 

2.17 

0.0004 

8.0 

2.17 

0.0004 

♦Relative  dielectric  constants  obtained  by  interpolation. 
♦♦Unable  to  measure. 

♦♦♦From  A.  R.  von  Hippel. 
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Table  D2 


Relative  Dielectric  Constant  and  Loss  Tangent  Data  for  Methanol* 

t 

I 


Frequency 

(GHz)  e 

r 


_ tan  6 _ 

Published**  Measured  Corrected 


3.0 

23.9 

0.640 

3.5 

4.0 

24.0 

0.218 

0.234 

4.5 

24.0 

0.199 

0.215 

5.0 

24.0 

0.173 

0.189 

I  Table  D3 

Relative  Dielectric  Constant  and  Loss  Tangent  Data  for  Distilled  Water* 


Frequency 

(GHz) 

e 

r 

tan  6 

Published** 

Measured 

Corrected 

2.5 

76.8 

0.136 

0.152 

3.0 

76.7 

0.157 

0.141 

0.157 

3.5 

73.6 

0.147 

0.163 

4.0 

70.9 

0.126 

0.142 

4.5 

68.5 

0.111 

0.127 

5.0 

66.4 

0.087 

0.104 

5.5 

64.5 

0.071 

0.087 

6.0 

62.7 

0.064 

0.080 

6.5 

61.1 

0.058 

0.074 

7.0 

59.6 

0.062 

0.078 

*ReLative  dielectric  constants  were  obtained  by  interpolation. 

**From  A.  R.  von  Hippel. 
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APPENDIX  E: 


ATTENUATION  OF  ELECTROMAGNETIC  WAVES  BY  A  HONEYCOMB  GRATING 


Data  showing  the  theoretical  attenuation  resulting  from  a  honeycomb 
grating  placed  in  a  pipe  are  presented.  The  amount  of  attenuation  depends  on 
the  grating  size  and  the  relative  dielectric  constant  of  the  fluid  in  the 
pipe.  A  problem  not  addressed  is  the  decrease  in  flow  rate  of  the  fluid  as  a 
result  of  the  grating. 

A  computer  program  was  developed  to  calculate  the  attenuation  of  an 
electromagnetic  wave  caused  by  a  grating.  The  incident  wave  is  taken  to  be 
normal  to  the  grating.  In  addition  to  the  attenuation  caused  by  the  waveguide 
below  cutoff  propagation,  the  attenuation  caused  by  scattering  at  the 
interfaces  is  computed.  Equation  Cl  in  Appendix  C  is  used  to  compute  the  loss 
from  waveguide  below  cutoff  propagation.  The  calculation  of  the  loss  due  to 
scattering  is  complex,  requiring  spectral-domain  analysis*^  and  will  not  be 
discussed  here. 

Figure  El  is  a  drawing  of  the  grating  used  in  the  analysis.  The  analysis 
was  done  for  four  fluids:  carbon  tetrachloride,  methanol,  distilled  water, 
and  air.  Tables  El  through  E4  list  the  attenuation  result  from  scattering  and 
waveguide  below  cutoff  propagation  for  various  frequencies  between  2.5  and 
7.6  GHz. 

Note  that  the  attenuation  from  scattering  is  minimal  and  can  be 
neglected.  This  suggests  a  simpler  method  for  analyzing  the  effects  of 
gratings.  Since  most  of  the  attenuation  provided  by  a  grating  is  caused  by 
waveguide  below  cutoff  propagation,  the  computer  programs  presented  in 
Appendix  C  can  be  used  to  compute  an  attenuation  factor.  Because  the  grating 
thickness  is  usually  much  less  than  a  wavelength,  the  single  mode  analysis 
used  by  programs  Pipe  1A  and  Pipe  IB  is  only  an  approximation.  Actual 
attenuation  could  be  slighly  higher  or  lower,  depending  on  the  amount  of 
energy  coupled  into  each  mode. 

In  the  simplest  case,  a  grating  can  be  viewed  as  a  collection  of 
waveguides.  It  is  important  to  note  that  in  this  approximation,  the  total 
attenuation  caused  by  the  grating  is  equal  to  the  attenuation  from  one  cell 
(waveguide) — not  to  the  sum  of  the  individual  cell  attenuations. 


^R.  Mittra,  R.  Hall,  and  C.  H.  Tsao,  "Spectral -Domain  Analysis  of  Circular 
Patch  Frequency  Selective  Surfaces  (FSS),"  IEEE  Trans.  Antenna  Propagation, 
in  press. 
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Figure  El.  Diagram  of  the  honeycomb  grating  used  in  calculating  the 
attenuations  in  Tables  El,  E2,  E3,  and  E4. 


Table  El 

Calculated  Attenuation  Factors  Resulting  Prom  a  Honeycomb  Grating: 

Carbon  Tetrachloride 


Frequency 

(GHz) 


Attenuation  from 
scattering 


Attenuation  from 
waveguide  below 
cutoff  propagation 


106.1 

106.0 

105.8 
105.6 

105.4 

105.1 

104.8 

104.5 

104.2 

103.8 
103.4 


56 


Table  E2 


B 

1 

B 

» 

9, 

2 


Calculated  Attenuation 

Frequency 

(GHz) 

Factors  Resulting  From 

Attenuation  from 
scattering 

a  Honeycomb  Grating:  Methanol 

Attenuation  from 
waveguide  below 
cutoff  propagation 

4.0 

.044 

96.4 

4.5 

.055 

93.6 

5.0 

.068 

Table  E3 

90.3 

Calculated  Attenuation  Factors  Resulting 

Distilled  Water 

From  a  Honeycomb  Grating: 

Attenuation  from 

Frequency 

Attenuation  from 

waveguide  below 

(GHz) 

scattering 

cutoff  propagation 

2.5 

.015 

93.6 

3.0 

.024 

87.5 

3.5 

.033 

79.9 

4.0 

.044 

70.3 

4.5 

.055 

57.9 

5.0 

.068 

40.6 

5.5 

.083 

16.1 

6.0 

.099 

11.1 

6.5 

.116 

6.4 

7.0 

.135 

6.4 

7.5 

.155 

5.2 

7.6 

.159 

Table  E4 

5.0 

Calculated  Attenuation  Factors  Resulting  From  a  Honeycomb  Grating:  Air 

Attenuation  from 

Frequency 

Attenuation  from 

waveguide  below 

(GHz) 

scattering 

-  cutoff  propagation 

2.5 

.015 

106.3 

3.0 

.024 

106.2 

3.5 

.033 

106.2 

4.0 

.044 

106.1 

4.5 

.055 

106.0 

5.0 

.068 

105.8 

5.5 

.083 

105.7 

6.0 

.099 

105.6 

6.5 

.116 

105.4 

7.0 

.135 

105.2 

7.5 

.155 

105.1 

_• 
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APPENDIX  F: 


TRANSMISSION  OF  ELECTROMAGNETIC  WAVES  THROUGH  MULTIPLE  DIELECTRICS 


Consider  the  problem  of  wave  propagation  through  a  multilayered 
dielectric  region  (Figure  FI).  Each  layer  has  relative  dielectric  constant 
(possibly  complex)  e  =  £.,  i  -  1,  2,  ...,  N  and  lies  between  z  =  z^_^  and  z^, 
where  Zq  =  -«  and  z^  =  +<=.  A  TEM  wave,  propagating  in  the  +z  direction,  is 
incident  from  the  left.  In  the  ith  region,  the  electric  and  magnetic  fields 
(E^  and  H^)  are  given  by 


where 


-jk.z  +  jk.z 

E.  *  T.  (e  1  +R.e  l) 

1  1  -jk.z  1  +  jk.z 

H.  =  T.Y.  (e  1  -  R.e  1 

ill  l 


) 


[Eq  FU 
[Eq  F2 ] 


Y.  =  /FT  Y  [Eq  F3 ] 

110 

k.  =  /FT  k  [Eq  F4] 

110  n 


Yq  and  kQ  are  the  free-space  wave  impedance  and  wave  number,  respectively.  T^ 
is  the  transmission  coefficient  in  the  ith  region.  R^  is  the  ith  reflection 
coefficient.  Since  regions  1  and  N  are  unbounded,  T^  =  1  and  RN  =  0. 


Forcing  the  continuity  of  E  and  H  at  the  interface  between  the  ith  and 
the  (i  +  l)th  layer,  the  following  relationship  results: 


-jk. z. 

J  l  l 


Y.  e 

l 


-jk.z. 
J  i  i 


w 


If  we  then  define: 


~Jki+lzi 


-jk. ^.z. 
J  l  +  l  l 


+ jk. z . 
11 


♦jk. z . 

li 


-Yi e 


+jk.  ,z. 
J  l  +  l  i 


+  jk 


Yi  +  le 


i  +  lZi 


T. 

l 

T.  R. 
l  i 


T. 


i  +  1 
Ti+1  Ri+1 


i  =  1,  2,  ...,  N-l 


[Eq  F5 ] 


X. 


[Eq  F6 ] 


A.  . 
l»J 


- jk. z . 
,  1  J 


Y.e 

i 


- jk. z . 
i  J 


+ jk. z . 
,  1  J 


-Y.e 

l 


+ jk. z . 
i  J 


[Eq  F7  ] 
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Figure  FI.  TEM  wave  propagation  through  a  muLtilayered  dielectric  region. 


Equation  F4  can  be  rewritten  as: 


Solving  for  X^: 


Ai,i  Xi  Ai+l,i  Xi  +  1  1  ~  1»  2»  ***»  N_1 


(Ai*l,i>  Ai,i  Xi  "  Xi+] 


[Eq  F8 ] 


[Eq  F9] 


Equation  F9  relates  the  reflection  and  transmission  coefficients  in  the  ith 
region  to  the  reflection  and  transmission  coefficients  in  the  (i+l)th  region 
by  a  matrix  equation.  These  matrices  can  be  cascaded  to  relate  X^  to  X^. 


XH  =  ^  (A.  .  .  )  A.  .  X. 

N  i=l  i 


[Eq  F10 ] 


The  bracketed  quantity  in  Equation  F10  is  a  2x2  matrix,  depending  only  on 
known  quantities.  As  Tj=l  and  R^=0  for  the  dielectric  slab  measurement, 
Equation  F10  consists  of  two  equations  and  two  unknowns,  R^  and  T^,  the 
overall  reflection  and  transmission  coefficients.  Equation  F10  can  be  solved 
for  Rj  and  T^.  The  matrix  product  in  Equation  F10  can  be  carried  out 
analytically,  although  for  large  N  values,  this  is  difficult.  It  is  often 
more  convenient  to  calculate  the  matrix  product  numerically. 

When  one  or  more  of  the  regions  is  particularly  lossy  and  is  more  than 
several  wavelengths  long,  certain  numerical  problems  can  arise  in  computing 
Equation  F10.  In  this  case,  the  A^  •  matrices  in  the  ith  layer  consist  of  both 
very  large  and  very  small  numbers.  To  evaluate  the  matrix  product  in  Equation 
F10,  these  numbers  must  be  multiplied  and  added.  This  can  yield  inaccurate 
results  when  using  finite  precision  arithmetic.  This  problem  can  be  avoided 
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Th«  Affect  of  fluids  on  waveguides  below  cutoff  penetrations  as  related  to 
electromagnetic  shielding  effectiveness  /  by  Michael  K.  Mclnemey  ... 
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tory,  1984. 
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1.  Shielding  (electricity).  2.  Wave  guides.  ).  Fluids.  L.  Mclnemey 
Michael  K.  11.  Ray,  Scott,  III.  McCormack,  Raymond  (i.  IV.  Castillo,  Stove 
V.  Mlctra,  Raj.  VI.  Series;  Technical  report  (Construction  Engineering 
Research  Laboratory)  M-354. 


